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Abstract 

We develop a method for detecting the inhomogeneity of the pion-emitting sources produced in 
ultra-relativistic heavy ion collisions, through event-by-event two-pion Bose-Einstein correlations. 
The root-mean-square of the error-inverse-weighted fluctuations between the two-pion correlation 
functions of single and mixed events are useful observables for the detection. By investigating 
the root-mean-square of the weighted fluctuations for different impact parameter regions people 
may hopefully determine the inhomogeneity of the particle-emitting in the coming Large Hadron 
Collider (LHC) heavy ion experiments. 
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Two-pion Hanbury-Brown-Twiss (HBT) interferometry is a useful tool for probing the 
space-time structure of the particle-emitting sources produced in high energy heavy ion 
collisions l|, |2|, y, |4j . Because of the limitation of data statistics, the usual HBT investigations 
are performed for mixed events and the HBT radii are obtained by fitting the two-pion HBT 
correlation functions from the mixed events to Gaussian parametrized formulas. 

On event-by-event basis, the density distribution of the source may be not a Gaussian 
distribution. An inhomogeneous particle-emitting source on event-by-event basis may be a 
more general case because of the fluctuating initial matter distribution in high energy heavy 
0, Q, B]. In Ref. jsj], a granular source model was used to explain the 



ion collisions 
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13|. The 



Relativistic Heavy Ion Collider (RHIC) HBT results, R ou t/R S ide 
main idea of the explanation is that the evolution time for the granular sources is short and 
the short evolution time, which can not be averaged out after event-mixing, may lead to 
the HBT results R out ~ -R S ide 



, [15I . Recent source imaging researches for the collisions 



of 



200 GeV Au+Au indicate that the pion-emitting source with the selections 



40 < centrality < 90% and 0.20 < < 0.36 GeV/c is far from a Gaussian distribution 
jltl [17I . Q • Although the long tail of the two-pion source function at large separation r is 



believed mainly the contribution of long-lived resonances 



16 



171 ]. the enhancement of the 
source function at small r may possibly arise from the source inhomogeneity fis| . 

For inhomogeneous particle-emitting sources, the single-event two-pion HBT correlation 
: unctions may exhibit fluctuations relative to the HBT correlation functions of mixed events 
19] . I20I ] . Detecting and investigating this event-by-event fluctuations is important for finally 
determining the source inhomogeneity and understanding the initial conditions and evolution 
of the system in high energy heavy ion collisions. 

Hydrodynamics may provide a direct link between the early state of the system and 
final observables and has been extensively used in high energy heavy ion collisions. In 
hydrodynamical calculations the system evolution is determined by the initial conditions 
and the equation of state (EOS) of the system. Smoothed Particle Hydrodynamics (SPH) 
is a suitable candidate that can be used to treat the system evolution with large fluctuating 

t has been used in 




25|. In 



initial conditions for investigating event-by-event attributes [7^ 
high energy heavy ion collisions for a wide range of problems 
the present letter we use SPH to describe the system evolution. The system initial states 
are given by the NEXUS event generator [6] at tq = 1 fm/c for y/s NN = 200 GeV Au+Au 
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FIG. 1: Transverse distributions of energy density of the systems with different impact parameters 
and at t = 1 fm/c [(a)-(c)] and t = 5 fm/c [(a')-(c')]. The unit of energy density is GeV/fm 3 . 

collisions at RHIC and tq = 0.5 fm/c for y/s NN = 5500 GeV Pb+Pb collisions at LHC The 
EOS is obtained with the entropy density suggested by QCD lattice results 
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In the 
gluons 



OS, the QGP phase is considered as an ideal gas of massless quarks (u, d, s) and 

the resonances with mass below 2.5 GeV/c 2 , 



22 



231 ] . The hadronic gas is composed o 



23j . The transition temperature between 



where volume correction is taken into account [22 
the QGP and hadronic phases is taken to be T c = 160 MeV, and the width of the transition 
is taken to be 0.1T r 



The coordinates used in SPH are r = \/t 2 — z 2 , x, y, and r\ = (1/2) ln[(t + z)/{t — z) 



22] . They are convenient for describing the system with rapid longitudinal expansion. 



However, in order to investigate the whole space-time structure of the system an nonlocal 
coordinate frame is needed and we work in the center-of-mass frame of the system. Figure 
1(a), (b), and (c) show the pictures of the transverse distributions of energy density for one 
y/s NN = 200 GeV Au+Au event at t — 1 fm/c and with impact parameter b = 0, 5, and 
10 fm, respectively. The pictures are taken for the spatial region (|x,?/| < 12 fm, \z\ < 0.83 
fm) and with an exposure of At = 0.3 fm/c. Fig. l(a'), (b'), and (c') are the corresponding 
pictures taken at t = 5 fm/c with the same At and for the spatial region (\x, y\ < 12 fm, 
\z\ < 1 fm). One can see that the systems are inhomogeneous in space and time. There are 
many "lumps" in the systems and the number of the lumps decreases with impact parameter 



increasing. 
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FIG. 2: Two-pion correlation functions for a sample of different single events (dashed lines) and 
mixed events (solid lines) with impact parameters b =10 fm (up panels) and b =5 fm (down panels). 



Assuming that final identical pions are emitted at the space-time configuration charac- 
terized by a freeze-out temperature Tf, we may generate the pion momenta according to 
Bose-Einstein distribution and construct the single-event and mixed-event two-pion corre- 
lation functions 0, Q. Figure 2 shows the two-pion correlation functions C(g si d e , (/out, <?ion g ) 
for the single and mixed events for ^ t /s^ = 200 GeV Au+Au with impact parameters b = 
10 fm (up panels) and b =5 fm (down panels). Here g S id e , q ut, and qi nT )g a re the components 
of "side", "out", and "long" of relative momentum of pion pair 29|, |30j. In each panel of 
Fig. 2, the dashed lines give the correlation functions for a sample of different single events 
and the solid line is for the mixed event obtained by averaging over 40 single events. In our 
calculations, the freeze-out temperature is taken to be 150 MeV. For each single event the 
total number of generated pion pairs in the relative momentum region (g S id e , <?out, Qiong < 200 
MeV/c) is N nn = 10 6 and the numbers of the pion pairs in the relative momentum regions 
(<?i < 200MeV/c; qj, < 30MeV/c) are about 2.7%N n7T , where i, j, and k denote "side", 
"out", and "long". 

From Fig. 2 it can be seen that the correlation functions for the single events exhibit 
fluctuations relative to those for the mixed events. These fluctuations are larger for bigger 
impact parameter. It is because that the number of the lumps in the system decreases 
with impact parameter and the fluctuations are larger for the source with smaller number of 
lumps 19|, |20j. It also can be seen that in the longitudinal direction the correlation functions 
exhibit oscillations which can not be smoothed out by event mixing. It is because that there 
are two sub-sources moving forward and backward the beam direction. This oscillations of 
the mixed-event correlation functions will not appear if we apply an additional cut for the 
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FIG. 3: The distributions dN/f for 40 events with FIC and SIC. b = 5 fm. 

initial rapidity of the "smoothed particles", t]q > or rjo < 0. 

We have seen that the two-pion correlation functions of the single events exhibit event- 
by-event fluctuations. However, in the usual mixed-event HBT measurements, these fluctu- 
ations are smoothed out. In order to observe the event-by-event fluctuations, we investigate 
the distribution dN/df of the fluctuations between the correlation 
mixed events, |C S (<&) — C m ((Zi)|, with their error-inverses as weights 



'unctions of single and 

3, 



Tm A\C s ( qi )-C m ( qi )\- {) 

In calculations we take the width of the relative momentum bin to be 10 MeV and use 
the bins in the region 20 < qi < 200 MeV. The up panels of Fig. 3 show the distributions 
of / in the "side", "out", and "long" directions, obtained from 40 simulated y/s NN = 200 
GeV Au+Au events. The impact parameter for these events is b = 5 fm and the number of 
correlated pion pairs for each of these events is = 10 7 . The solid lines are the results for 
the events with the fluctuating initial conditions (FIC) generated by NEXUS. Because in the 
last analysis the fluctuations of the single-event correlation functions are from the FIC, for 
comparison we also investigate the distribution dN / df for the events with the smooth initial 
conditions (SIC) obtained by averaging over 30 random NEXUS events jl, [si, 0, 24] with 
the same b as for the FIC. It can be seen that the / distributions for FIC are much wider 
than the corresponding results for SIC. The down panels of Fig. 3 show the distributions 
of / for the 40 simulated events and N nn = 5 x 10 6 . It can bee seen that the widths of the 
distributions for FIC decrease with N nn decreasing and for N nn = 5 x 10 6 the distributions 
for FIC are still wider than those for SIC. 

In experiments the number of correlated pion pairs in one event, N nn , is limited. It is 
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FIG. 4: The distributions dN/f for 40 events with FIC and SIC, where / are calculated for variables 
gtmns and q. 

related to the energy y/s NN of the collisions. For a finite N nn we have to reduce variable 
numbers in analysis although it will lose some details. In Fig. 4 we show the distributions of 
/ for the variables of transverse relative momentum gtrans and relative momentum q of the 
pion pairs for the 40 simulated events with b = 5 fm. One can see that for N n7T = 5 x 10 6 , 
the distributions for FIC are much wider than those for SIC both for gtrans arid q. Even for 
N n7T = 5 x 10 5 , the widths for FIC are visibly larger than those for SIC. In order to examine 
the distributions quantitatively, we calculate the root-mean-square (RMS) of /. Figure 5 
shows the RMS, / rms , as a function of N n7T for the 40 simulated events for y/s NN = 200 GeV 
Au+Au with b = 9 fm (the up panels) and b = 5 fm (the down panels), where the SIC are 
obtained by averaging over 100 NEXUS events. It can been seen that the values of / rms 
rapidly increase with N vn for FIC because the errors in Eq. (CQ) decrease with N„ n . For FIC 
the results for b = 9 fm are larger than the corresponding results for b = 5 fm because the 
differences \C s (qi) — C m (qi)\ in Eq. (TjQ) increase with b. For SIC the values of / rms are almost 
independent from N n7T . It is because that both the differences and their errors in Eq. (TjQ) 
decrease with N n7T in the SIC case. 

In Fig. 5 we only display the statistic errors which may decrease with the number of 
events. Moreover, in order to display the fluctuations of the single-event HBT correlation 
functions and the effect of FIC on the distributions of /, we used very large N nn in our 
calculations. At RHIC energy the event multiplicity of identical pions, M n , is about several 
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FIG. 5: The root-mean-square of / for variables q and cftrans as a function of N w 

2 

1.5 



3 _ (a) LHC, FIC 



co 

E 2 

1 




(c) RHIC, FIC 



(d) RHIC, SIC 




(b) LHC, SIC 
o— © fo=9±0.3fm 
A- -A fc=5±0.3fm 




1 

0.5 



1 10 50 1 10 50 

A^(xl0 5 ) A^(xl0 5 ) 



FIG. 6: The root-mean-square of / calculated with cjtrans for 40 simulated events at LHC and RHIC 
energies. 

hundreds for central collisions. The order of N nw is about 10 5 (~ Ml/2). However, at the 
higher energy of LHC, M w will be about two thousands and the order of N W7r will be 10 6 . In 
this case the large differences between the RMS of / for inhomogeneous and homogeneous 
sources provide a great opportunity to detect the source inhomogeneity. 

Another problem in experimental data analysis is that the impact parameter b is hardly 
to be held at a fixed value and usually be limited in a region. In such a case the large / rms 
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values may arise from the source inhomogeneity as well as the variation of b in the region. 
So one should also consider the effect of the b variation in inhomogeneity detections. 

Figure 6(a) and (b) show the RMS of / calculated with q tTans for the 40 simulated events 
for ^/s NN = 5500 GeV Pb+Pb with FIC and SIC, respectively. The impact parameters 
for the FIC events are selected randomly from the regions b = 9.0 ± 0.3 (circle symbol) 
and b = 5.0 ± 0.3 (triangle symbol), respectively. Whereas for each SIC event the SIC is 
constructed by averaging over 30 NEXUS events with the same impact parameter selected 
randomly from the b region as for the FIC case. Figure 6(c) and (d) give the corresponding 
results for ^/s NN = 200 GeV Au+Au collisions. Our calculations indicate that the lump 
numbers in the sources at the LHC energy are larger than those at the RHIC energy. So the 
RMS values for the LHC energy are smaller than those for the RHIC energy. From Fig. 6 
one can see that for FIC case (inhomogeneous sources) the RMS for the larger b are higher 
than those for the smaller b. It is because that the number of lumps in the inhomogeneous 
source decreases with b (see Fig. 1). However, for the SIC case the RMS for the two impact 
parameter regions are almost the same for fixed N nn . In this case the increases of the RMS 
values with N nn are due to the fluctuations of the correlation functions arising from the 
variation of b in their regions 5b. Our calculations indicate that these increases will become 
flat when 5b — > 0. So, one may determine the source inhomogeneity by analyzing the RMS 
values for different b regions. For inhomogeneous sources their differences are bigger for 
larger N wn . Otherwise, they are almost the same. 

In summary, on event-by-event basis the initial density distribution of matter in high 
energy heavy ion collisions is highly fluctuating. The fluctuating initial conditions lead 
to event-by-event inhomogeneous particle-emitting sources. In this letter we developed a 
method for detecting the source inhomogeneity through event-by-event two-pion correla- 
tions. We find that the RMS of the error-inverse-weighted fluctuations / are useful observ- 
ables for detecting the inhomogeneity of the sources. The high identical pion multiplicity in 
the coming LHC heavy ion collisions provides a great opportunity to do the detections. By 
investigating the RMS values of / for different impact parameter regions people may hope- 
fully determine the inhomogeneity of the particle-emitting produced in the coming LHC 
experiments. 
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